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We define the pion-photon Transition Distribution Amplitudes (TDA) in a field theoretic formal- 
ism from a covariant Bethe-Salpeter approach for the determination of the bound state. We apply 
our formalism to the Nambu - Jona Lasinio model, as a realistic theory of the pion. The obtained 
vector and axial TDAs satisfy all features required by general considerations. In particular, sum 
rules and polynomiality condition are explicitly verified. We have numerically proved that the odd 
coefficients in the polynomiality expansion of the vector TDA vanish in the chiral limit. The role of 
I PCAC and the presence of a pion pole are explicitly shown. 
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Q\ • I. INTRODUCTION 

(N ■ 

Hard reactions provide important information for unveiling the structure of hadrons. The large virtuality, Q^, 
• involved in the processes allows the factorization of the hard (perturbative) and soft (non-perturbative) contributions 
Oh' in their amplitudes. Therefore these reactions are receiving great attention by the hadronic physics community. In 
^ . the past, only total cross sections of inclusive processes or longitudinal asymmetries, that have simple parton model 
' interpretations, were studied. The basic theoretical ingredients to be understood are diagonal Parton Distribution 
^ ■ Functions (PDF) Q, governing the deep inclusive processes. In the recent years a new variety of processes, like the 
Deeply Virtual Compton Scattering, has been considered. These processes are governed by the Generalized Parton 
Distributions (GPD) 0, i, i, [|. From a theoretical point of view, PDFs are related to diagonal matrix elements 
j>! : of a bilocal operator from the initial hadron state to the same final hadron state (the same particle with the same 
\Q momentum). The GPDs are related to the matrix elements of the same bilocal operator as in the previous case, where 
\C the initial and final hadron are the same particle, but have different momentum. The GPDs describe non-forward 
' matrix elements of light-cone operators and therefore measure the response of the internal structure of the hadrons 
, to the probes. 

t — ■ The generalization of parton distributions to the case where the initial and final states correspond to different 
] particles has recently been proposed in @, Such distributions are called Transition Distribution Amplitudes 
■ (TDA) since they have been introduced through hadron-photon transitions. In particular, the easiest case is to 
consider pion-photon TDA, governing processes like tt+tt" 7*7 or 7*7r+ — > 77r+ in the kinematical regime where 
the virtual photon is highly virtual but with small momentum transfer. 

The aim of this work is to calculate the pion-photon TDAs in a field theoretic scheme treating the pion as a bound 
state in a fully covariant manner using the Bethe-Salpeter equation. In this way we preserve all invariances of the 
problem. In order to perform a numerical study we will use the Nambu - Jona Lasinio (NJL) model to describe the 
pion structure. The NJL model is the most realistic model for the pion based on a local quantum field theory built 
with quarks. It respects the realizations of chiral symmetry and gives a good description of the low energy physics of 
the pion [1]. 

The NJL model is a non-renormalizable field theory and therefore a cut-off procedure has to be defined. We have 
chosen the Pauli-Villars regularization procedure because it respects all the symmetries of the problem. The NJL 
model together with its regularization procedure is regarded as an effective theory of QCD. Moreover, it has been 
used to tune coefficients of Chiral Perturbation Theory ^9'] . 

The NJL model has been used to describe the soft (non perturbative) part of the deep processes, while for the hard 
part conventional perturbative QCD must be used. It has been applied to the study of pion PDF [l^, [HI, [13 and 
to the pion GPD [ij]. In the chiral limit, its quark valence distribution is as simple as q{x) = (x) 9{1 — x). Once 
evolution is taken into account, a good agreement is reached between the calculated PDF and the experimental one 
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FIG. 1: Pion pole contribution between the axial current (represented by a cross) and the photon-external pion vertex associated 
to the last contribution of equation ([2]). 



[lO| . A more elaborated study of pion PDF is done in Ref. [l4| using non local lagrangians [l5| , which confirms that 
the result obtained in the NJL model for the PDF is a good approximation. 

As GPDs, TDAs must satisfy sum rules. The vector and axial TDA are connected to the vector and axial form 
factors, Fy and Fa, appearing in the tt"*" — > ^e'^v process. In order to have a proper understanding of the axial 
TDA we will make the difference between two contributions to the axial current coming from the analysis of the 
amplitude for the pion radiative decay. The first one is originated from the internal structure of the hadron, in our 
case a pion. This contribution is the proper axial TDA. A second contribution is present and it can be understood as 
a manifestation of PCAC. The axial current can be coupled to a pion. This pion is a virtual one and its contribution 
will be present independently of the external hadron. 

The polynomiality condition is also satisfied. We observe that, in the polynomial expansion of the moments of the 
vector TDA, the odd powers of ^ are chirally suppressed and that they vanish in the chiral limit. 

Previous studies of the axial and vector pion-photon TDA have been done using different quark models [l^, [13] . 
Both works parametrize the TDAs by means of Double Distributions. 

This paper is organized as follows. In section |TT] we establish the connection between the TDA and the vector and 
axial pion form factors, Fy and F^- In section Hill we define our approach for the TDA and we calculate them in the 
NJL model. In section IIVI we study the sum rules and the polynomiality condition of the TDAs. In section |V] we 
discuss our results and we finally give our conclusions in section IVll 



II. THE PION-PHOTON TRANSITION DISTRIBUTION AMPLITUDES 



The pion-photon TDAs are connected, through sum rules, to the vector and axial-vector pion form factors, Fy and 
Fa- Before giving a proper definition of the TDAs let us recall the definition of these form factors. They appear in 
the vector and axial vector hadronic currents contributing to the decay amplitude of the process 7r+ "/e^v. The 
precise definitions of these currents are [H, [l^ : 

{l{p')\q{0)l^.r-q{0)\7T{p))^~tee-'e^,,,p'Pp'^^^ , (1) 



(7 (P')l 9 (0)7^75^-9(0) kb)) ^ee-^ {p'^P. - g^.p' .p) 



+ eeM(y-p)„p.^4^-A/2/..9,,J , (2) 



with /tt — 93MeV, e^^^s _ ^ a^cj — (j^ _ ^7-2) j2. All the structure of the decaying pion is included in the form 
factors Fy and Fa- We observe that the vector current only contains a Lorentz structure associated with the Fy 
form factor. The axial current is composed of two terms. The first one, defining Fa, gives the structure of the pion. 
The second one corresponds to the axial current for a point-like pion. It has two different contributions. The first 
one corresponds to a point-like coupling between the incoming pion, the outcoming photon and a virtual pion which 
is coupled to the axial current. It is depicted in the diagram of Fig. [T]and can be seen as a result of PCAC, because 
the axial current must be coupled to the pion. It isolates the pion pole contribution of the axial current in a model 
independent way. The second contribution of this term, proportional to j-Kg^iv, corresponds to a pion-photon-axial 
current contact term. With these definitions, all the structure of the pion remains in the form factor Fa- 

Let us go now to TDAs. For their definition we introduce the light-cone coordinates — (w" ± w'^) j \/2 and the 
transverse components = for any four-vector - We define P = (p + p') /2 and the momentum transfer, 

A = p' — p, therefore = — t/4 and t = A^. The skewness variable describes the loss of plus momentum 
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of the incident pion, i.e. ^ = {p — p') /ZP"*", and its value ranges between t/ (2rn^ — t) < ^ < 1. Actually there 
is no symmetry relating the distributions for negative and positive ^ which could have constrained the values of the 
skewness variable to be positive, like for GPDs. The vector and axial TDAs are the Fourier transform of the matrix 
element of the bilocal currents, g (— z/2) 7^ [75] g (z/2), separated by a light-like distance. Then, they are directly 
related to the currents defined in Eqs. ((T][2l) through the sum rules: 



=0 



qm^{p)) ■ (3) 



With this connection we can introduce the leading twist decomposition of the bilocal currents. For that we introduce 
the light-front vectors p^" = P+ (1, 0, 0, 1) /a/2 and nf" = (1, 0, 0, -1) / (\/2P+) . The explicit expressions for the pion 
and photon momenta in terms of their light-cone components are given by Eqs. (|Aip and (|A2p . Then we have 



P+ 



ee e 



pP 



V2U 



(4) 



l^e-^- (7 (.01 (-1) ^a.-. (0 k ip))l__^^__^ - • A^) 



V2U 



(5) 



where e {£,) is equal to 1 for ^ > 0, and equal to —1 for ^ < 0. Here V (x, ^, t) and A {x, ^, t) are respectively the vector 
and axial TDAs. They are defined as dimensionless quantities. From the condition ^ we observe that they obey the 
following sum rules: 



' dx V^* ix,^,t) = ^^Fv{t) , 
1 rrij, 



L 

pi 

j dx A'^'' (x,^,t) 



V2U 



FA{t) , 



(6) 



(7) 



which were firstly introduced in Refs. [1,0|- 

In the second term of Eq. (O, we have introduced the Pion Distribution Amplitude (PDA) </> (x) . By definition the 
PDA is 



27r 



= ■^iV2U(Pix) 



(8) 



The PDA vanishes outside the region x £ [0, 1] and satisfies the normalization condition 



dx (j) {x) = 1 



(9) 



This second term has been introduced in order to isolate the pion pole contribution of the axial current in a model 
independent way, as we have done in Eq. ^ for the 7r+ —t ^e^u process. Therefore, all the structure of the pion 
remains in the TDA A {x, t). It can be seen as a result of PCAC, because the axial current must be coupled to the 
pion. Therefore, this term is not a peculiarity of the pion-photon TDAs. A similar pion term will be present in the 
Lorentz decomposition in terms of distribution amplitudes of the axial current for any pair of external particles. A 
pion exchange contribution has already been analyzed in [2O, [2l| for the axial helicity-flip GPD and, in [3| : a similar 
structure for the axial current has been obtained using different arguments^. This term we have represented in Fig. [2] 
is only non-vanishing in the ERBL region, i.e. the x € [—£,, S,] region. The kinematics of this region allow the emission 
or absorption of a pion from the initial state, which is described through the PDA. And it can be seen from Fig. [2] 
that positive values of 1^ corresponds to an outcoming virtual pion, whereas negative values of ^ describe an incoming 
virtual pion. The latter is related to the matrix element (tt {p) \ qfi^^r^q |0) , instead of the one present in Eq. ([H]), 
what gives rise to the minus sign included in e (^) . 



^ In order to make this connection it must be realized that there is a \/2 between our definition of fn and the one used in [T^ . and that 
e-^.A^ = (1 -0 (e-A) 
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FIG. 2: Pion pole contribution to the axial bilocal current corresponding to the last term of equation ([5]). 




FIG. 3: Diagrams contributing to the TDA. We have depicted diagrams in which a quark u is change into a quark d by the 
bilocal current. There are similar diagrams in which the antiquark d is changed into a u 



III. A FIELD THEORETIC APPROACH TO THE PION-PHOTON TDA. 

In Ref. [l^ we have defined a method of calculation for the pion GPD in a field theoretical scheme, treating the 
pion as a bound state of quarks and antiquarks in a fully covariant manner using the Bethe-Salpeter equation. We 
apply here the same method for evaluating the pion-photon TDAs. This method has enormous advantages because 
it preserves all the physical invariances of the problem. Therefore, any property as sum rules or polynomiality is 
preserved. 

As usual, we consider that the process is dominated by the hand-bag diagram. Each TDA has two related contribu- 
tions, depending on which quark (u or d) of the pion is scattered off by the deep virtual photon. In Fig. [3]we depicted 
the diagrams in which the photon scatters off the u-quark. We observe that there are two kinds of contributing 
diagrams. In the first one the d-antiquark appears as the intermediate state, while in the second the bi-local current 
couples to a quark-antiquark pair coupled in the pion channel. The latter is present only for the axial current and 
includes the pion pole contribution. 

The details of the method of calculation are given in Ref. [13] ■ In the present case we obtain, from the first kind 
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of diagram of Fig. [3l the following contributions 
dz 



^e-^^^- (7 ip', e) \q (-|) 7^ [75] r-q (|) | 7r+ (p) 



^ ( J + T3 ) 7^^ ip' - k) r [75] T- iS ip - k) {k,p) iS i-k) 



Tr 



- 5 ( xP+ - - (2fc+ -p+ - p'+] 



Tr 



(10) 



where S (p) is the Feynman propagator of the quark and (f)^ {k,p) is the Bethe-Salpeter amplitude for the pion. Here 
Tr() represents the trace over spinor, color and flavor indices. The first contribution in Eq. pil)) is the one depicted 
in the first diagram of Fig. [3l The second contribution corresponds to a similar diagram but changing quarks u and 
d. In the NJL model, (j)'^^ {k,p) is as simple as 

(/.^^(fc,p) = -5^,,75V^r+ , (11) 

where g^qq is the pion-quark coupling constant defined in Eq.f IA5) ). 

The vector TDA has contribution only from this first kind of diagrams. We can express V {x,S,,t) as the sum of 
the active u-quark and the active d-quark distributions. The first contribution will be proportional to the d's charge, 
and the second contribution to the w's charge. Therefore, we can write 



Isospin relates these two contributions, v'^^- (x,^,i) — v^^^ {—x,S,,t) . A direct calculation gives 

(X, i) = 8 Nc U gTrqq m (x, ^, t) , 

with 



d^k 



(/e^ — rn^ + ie) (p' — k) — m"^ + ie [p ~ k) — m? + ie 



(12) 



(13) 



(14) 



In this integral, we first perform the integration over k . The pole structure of the integrand fixes two non- vanishing 
contributions to v^^j^ {x,^,t), the first one in the region ^ < a; < 1, corresponding to the quark contribution, and 
the second in the region — ^ < x < corresponding to a quark- antiquark contribution. Given the relation ([T2|) . the 
support of the entire vectorial TDA, (x,^,i), is therefore x e [—1, 1]. The analytical expression for (fT4| is given 
by Eq. (jAll|l . For the tt", the contributions of the first type of diagram (Fig. [3|) can be related to w^^^ {x, ^, t) 



Vf = ^ (yll, {x, 1 1) + vll, i-x, t)) , 



V2 

vf = ^ (x, 1 1) + vll, i^x, t)) . 



(15) 



Turning our attention to the axial TDA, we find a new contribution arising from the second diagram of Fig. [31 
This second contribution comes from the re-scattering of a qq pair in the pion channel. Therefore 



Zn \ 



2+=Z^=0 



2»g 



(16) 
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where i is an isospin index, 

M' = -e 

-Tr 



-Tr 



1 /I 



6- {k',p) iS {-k')-l-+T3] Y ip' - k') Z75r' iS {p - k') 



{k',p) iS [k' - p) i^^T^ iS {k' - p') \(^-+ T3) 7" iS [k') I , 



(17) 



j -0^3 (^xP+ - i {p'+ +P+- 2k+)^ Tr [zS (p' - k) t^Tst" iS {p - k) i^^t^] , 



and Hps is the pseudoscalar polarization, 



lips (A^) 



(18) 



(19) 



We can now evaluate the axial current in a straightforward way. Nevertheless, in order to extract the axial TDA we 
must subtract the pion pole contribution. We need for that the pion amplitude which, in the NJL model, is 



mg^gg Ne , [m^ - mix (1 - x)] 



[mf — m^x (1 — a;) 



< a; < 1 



(20) 



where and to^ are defined in the Appendix. Now, after a long but direct calculation, we obtain the expression 
for A (x, ^, i) . As the vector TDA, A {x, ^, i) can be expressed as a sum of the contributions coming from the active 
u-quark and the active J-quark. The first one will be proportional to the rf's charge and the second to the m's charge. 



A^^ {x, i,t) = - ^all^ {x, t) + ^a|Xa t) 



(21) 



Isospin relates these two contributions, dj^^ {x,^,t) = —au^a{—x,^,t) , where the minus sign is originated in the 

change in helicity produced by the 75 operator. The expression for a^^^ {x, ^, t) depends on the sign of In the 
1^1 < a; < 1 region we have 



{x, t) = -8 U g^gg m\{l+ ""Hl/^+j^, .f^^!^f ) ^3. {x, t) 



2^ml-t{l + 1-^ 



+ 



1 + e 1 



log- 



m 



m? (m^ — zm'^) 



(22) 



with the abbreviations 2; = (1 — a;) (^ + x) / (1 + ^) . And in the — |^| < a; < |^| region we have 
all, {X, ^,t) = -8N^ U g^gg m | (l + ""^^^^ - ^(1 + O^"^ ^ I^) ^' 



+ 



+ 



e(0 1 
1 - ^ 167r2 -f^ 



E 



1+e 



■log 



log 



(4e 



(4^^m^ — a; (m? - y ml) 
t) {4:^'^mf -xml 



^m^ - X 



t — ml [AS^'^mf — xt) [AS^'^m^ — x ml) 



(23) 



where a; = (^^ — x^) and y = 2; for ^ > and ?/ = for ^ < 0. 

The expressions for both the vector and axial TDAs in the chiral limit, i.e. m^- = 0, are well defined. In particular, 
for t = 0, we find the following simple expression for v'!^,{x,£,,t) 



vll,{x,^,0) = y/2U6F^\0) 



2|el(l + iei) 



(24) 



where Fy ^(0) is the chiral limit of the vector pion form factor at zero momentum transfer Fy -^(0) = 
lim„,^o-Fy^(0)/TO^ = 0.17 GeV-^ 
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n 


1 


2 


3 


4 




Ci.o(t) 


C2.0 (t) C2,l (t) 


Ca,o (t) C3.1 (t) C3,2 (t) 


C4,0 (t) Ci,l (t) C4,2 (t) C4.3 (t) 


nirr = 










t = 
t = -0.5 
t = -1.0 


22.6 
13.7 
10.4 


-22.6 0.0 
-16.3 0.0 
-13.4 0.0 


3.77 0.0 3.77 
3.00 0.0 2.43 
2.57 0.0 1.90 


-6.7 0.0 -6.7 0.0 
-5.7 0.0 -4.9 0.0 
-5.0 0.0 -4.0 0.0 


m.^ = 140 










t = 
t = -0.5 
t = -1.0 


22.7 
13.7 
10.3 


-22.9 0.44 
-16.4 0.25 
-13.4 0.19 


3.80 -0.10 3.67 
3.00 -0.06 2.43 
2.57 -0.05 1.87 


-6.8 0.19 -6.8 0.12 
-5.7 0.13 -4.9 0.08 
-5.0 0.11 -4.0 0.06 



TABLE I: Coefficients of tfie polynomial expansion for the vector TDA. The pion mass is expressed in MeV and t is expressed 
in GeV^. Notice that the coefficients have to be multiplied by 10~^. 



A similar expression is obtained for a'^^j^ix, ^, t) in the chiral limit for t = 



4^2(1 + 1^1) 



(l + l^l) 



(1-^(1-6 



(25) 



with '^(O) the axial form factor at i = in the chiral Umit ^(0) = Hm^^^o {0)/m^ = Fy ^(0). 



IV. SUM RULES AND POLYNOMIALITY 



The vector TDA of the 7r+ must obey the sum rule given in Eq. ([6]) with the expression of the vector pion form 
factor in the NJL model given by Eq. (|A7p . We have numerically recovered the sum rule for different t values. In 
particular we obtain the value ^^^(0) = 0.0242 for the vector form factor at t = 0, which is in agreement with the 
experimental value Fv{0) = 0.017 ± 0.008 given in [2^ . 

The 7r° distribution must satisfy the following sum rule Q 



1 

dx 

-1 



{QuV:"ix,C,t)-QdVfix,C,t))^V2UF^^,^it) . (26) 



A theoretical prediction of the tt" form factor value is given in [2^. In particular, at t = 0, the value F^.^»^(0) = 0.272 
GeV^^ is found. The neutral pion form factor is directly related to the 7r+ vector form factor so that the sum rule 
is satisfied. We obtain the value F7r^.^(0) = 0.244 GeV~^. In [23|, a dipole parametrization based on experimental 
data is proposed for the ^-dependence of ^.^-7*7 (t), obtaining for the dipole mass A = 0.77 GeV. We have found that, 
for small values of t, the NJL neutral pion form factor can be parametrized in a dipole form with A — 0.81 GeV. 

The axial TDA obeys the sum rules given by Eq. ([7]) with the axial form factor given in the NJL model by 
Eq. (|A9[) . This sum rule is satisfied for different t values. The numerical results also coincide. In particular we obtain 
FJ+(0) = 0.0239 for the axial form factor at t = 0, which is about twice the value FX^ (0) = 0.0115 ± 0.0005 given by 
the PDG [13]. 

We expect TDAs to obey the polynomiality condition. However no time reversal invariance enforces the polynomials 
to be even in the A-momenta like for GPDs. That means that the polynomials should be "complete" , i.e. they should 
include all powers in ^, 

/I n-l 
da;x"-iF(a;,e,t)-^ a^Wr • (27) 

However, in the chiral limit and for t — 0, we have analytically found that the odd powers in ^ go to zero for the 
polynomial expansion of the vector TDA. A study of the polynomiality in the limit given in Eq. (|24|) leads to the 
following analytical expression for the coefficients 2i (0 



ci^M^^U^K'^'io) (-i + 2(-i)"-i) (^-^ 



(28) 
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n 


1 


2 


3 




C'i,o (t) 




c!i,o{t) c^.iW c',,At) 


nin = 140 








t = 
t = -0.5 
t = -1.0 


22.4 
16.1 
13.2 


-3.77 4.00 
-2.97 2.57 
-2.53 1.97 


6.8 -4.9 2.3 
5.7 -3.5 1.7 
5.0 -2.8 1.3 



TABLE II: Coefficients of the polynomial expansion for the axial TDA. The pion mass is expressed in MeV and t is expressed 
in GeV^. Notice that the coefficients have to be multiphed by 10"'^. 



Notice they do not depend on i. We have numerically tested the polynomiality and obtained it. We observe that the 
coefficients for the odd powers in ^ are of one order of magnitude smaller than those for the even powers in ^. In 
particular we have numerically proved that, in the chiral limit, the polynomials only contain even powers in ^ for any 
value of t, 



(29) 



The coefficients in the chiral limit, Eq. (|29p. as well as the coefficients for m^r = 140 MeV, Eq. (|27p. we numerically 
obtained are given in Table IH 

The ^-dependence of the moments of the axial TDA A{x,^,t) also has a polynomial form. Those polynomials 
contain all the powers in ^, i.e. even and odd, 

/I n-l 
dxx--'A{x,^,t) = Y.^Ut)C ■ (30) 
■1 i=0 



An analytic study of the polynomiality in the limit given in Eq. (|25p confirms that all the powers in ^ have to be 
present. The analytic values for the coefficients in this specific limit are 

C'^%i0)^V2U2Ff^0) (l + 2(-l)«-i) ^"2* 



n{n + l){n + 2) 



C:':„«(0) = ^/2/. 2FfH0) (l + 2(-l)--) ^^^f '^^l . (31) 

which is in agreement with those numerically obtained. The coefficients of the polynomial expansions in the chiral 
limit are very close to those obtained for the physical values of the pion mass. In Table HIl the coefficients ^{t) are 
therefore shown only for ttItt = 140 MeV. The polynomiality property of the term containing the pion DA can also 
be studied. The i-dependence only comes from the pion pole. We can therefore write a general relation about the 
polynomiality property of the whole axial bilocal matrix element 

f^dxx--' I ^e"P"^-(^(P')|g-(-0 l,n^l,r-q[^)\niP)) = ^ e (e.A) ^ ^.(i) T ■ 

(32) 



V. DISCUSSION. 



In Figs. |3]and[ni the vector and axial TDAs are plotted in function of x for several values of ^ and t. In these figures, 
the it-quark contributes to the TDAs in the region of x going from — |^| to 1 and the d-antiquark going from —1 to |f | . 
Therefore, for x S [|^| ,1] {x E [—1, — |^|]) only w-valence quarks (d-valence antiquarks) are present (DGLAP regions). 
Besides, TDAs in the — |^| < x < |^| region (ERBL region) receive contributions from both type of quarks. For the 
vector TDA, we observe in figure |4] that the position of the maxima is given by the ^ value, separating explicitly the 
ERBL region from the DGLAP regions. The vector TDA is positive (negative) for negative (positive) values of x with 
the change of sign occurring around a; = 0. This change in the sign of the vector TDA is originated in the presence 
of the electric charge of the quarks in Eq. (fT2|) . 
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The process involved in the calculation of the TDAs allows negative values of the skewness variable. In the chiral 
limit, the skewness variable goes from ^ = — 1 to ^ = 1, for any value of t. In the chiral limit, the vector TDA for 
a negative value of f is equal to the vector TDA for |^|. This can be seen from the polynomiality expansion in this 
limit, Eq. (|29p. which has only even powers of 5. For the physical value of the pion mass, negative values of ^ are 
bounded by t/ (2m^ — t) < ^. For each allowed value of ^, we found that the numerical value of V (x, — |^| , t) is close 
to V (x, 1^1 , t) , due to the smallness of the coefficients of the odd powers of ^ in the polynomial expansion (|27|) . 

Analyzing the axial TDAs, plotted in Fig. [5l we observe two different behaviours depending on the sign of ^. For 
positive ^, the position of the minima is given by the value of the skewness variable while the position of the maxima 
is always x ~ in the ERBL region and x ~ ±(1 + ^)/2 in the DGLAP regions. For negative ^, the value of the axial 
TDA at X = ±^ is important and in some cases is a maximum and, in the ERBL region, A (x, ^, t) presents a minimum 
near x = 0. As we have previously shown, the axial TDA in the ERBL region receives contributions from two different 
diagrams, depicted in Fig. [3] In the second of these diagrams, a virtual quark-antiquark interacting pair in the pion 
channel appears. The pion pole, contained in this diagram, has been subtracted but the remaining non-resonant part 
contributes to the axial TDA. We observe that the latter contribution is the dominant one and produces the maxima 
around x = for positives f (Fig. (S]). Now, the axial TDA does not change the sign when we go from negative to 
positive values of x. In the axial TDA the change of sign originated in the presence of the electric charge of the quarks 
in equation (PTjl is compensated by the change of sign between quark and antiquark contributions generated by the 
75 operator present in the axial current. 

By comparing the plots for different values of the momentum transfer, it is observed that the amplitudes are lower 
for higher {—t) values, as it can be inferred from the decreasing of the form factors with {—t), connected to the TDAs 
through the sum rules. By increasing the {—t) value, not only the width and the curvature of the TDAs are changed 
but we also observe that higher values of ^ are preferred, i.e. the sign of the derivative of the collection of maxima 
changes passing from a zero momentum transfer to a non-zero one. 

Isospin relates the value of the vector and axial TDAs in the DGLAP regions, 

Vix,U) = -lv{-x,^,t), A{x,tt)^^A{-x,^,t), m<x<l , (33) 

being the factor 1/2 the ratio between the charge of the u and d quarks. We observe in Figs. 2] and [5] that our TDAs 
satisfy these relations. It must be realize that the relation ((33|) cannot be changed by evolution. 

Regarding the chiral limit, we observe that both the vector and axial TDAs do not significantly change going from 
a non-zero pion mass to the physical mass, except for the change in the lower bound of ^. 

Previous studies of pion-photon TDAs have already been done (igI. [iTj. In both works, double distributions have 
been used. Even if our model is different, a comparison of the results is still possible. The aim of the author of the 
first paper [lH is to provide some estimates of the vector and axial TDAs on the basis of the positivity bounds. The 
order of magnitude of the obtained amplitudes are similar to ours, but the former are constrained by the Fy and 
Fa form factors, through the sum rules. The sum rules are an input imposed in Ref. [l^ while it is a result in our 
calculation. The vector TDA obtained in this paper has some similarities with ours. Nevertheless, it does not satisfy 
the isospin relation (|33p due to the different choice of the M-quark and d-antiquark distributions, the first one related 
to the pion and the second to the photon, used in the saturation of the positivity bounds. We have also studied the 
positivity bounds for GPDs and noticed that it is actually an upper bound that is sometimes very higher than the 
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FIG. 5: The axial TDA in the case = 140 MeV and, respectively, t = and t = -0.5 GeV^. 




FIG. 6: Contributions to the axial TDA for both positive — 0.25, plain line) and negative = —0.5, dashed line) values 
of the skewness variable and for = 140 MeV and t — —0.5 GeV^. In each case, and in the ERBL region, the contribution 
coming from the first diagram of Fig. [3] is represented by the dashed-dotted lines and the non-resonant part of the second 
diagram of Fig. [3] is represented by the dotted lines. 



value of the GPD itself. The vector TDA obtained in [16] is rather peaked at a; = whereas we have a smoother 
behaviour. 

A detailed comparison with the results of Ref. [l3| is not easy, due to the choice of the asymmetric notation. 
Nevertheless the vector TDA seems to be consistent with our results and some limits, in particular when — and 
t = 0, can be recovered going from one notation to the other. 

Regarding our axial TDA, it differs from the two previous calculations [IB, [13 "^^^ effect of the non-resonant 

part of the second diagram of Fig. [3l This contribution, corresponding to the last term of Eq. ([23]) . is proportional 
to (t — TO^) but with zero value for the residue. The presence of this term is crucial in order to obtain the axial 
form factor using the sum rule. In fact, its contribution in Eq. (jA9p can be easily recognized. Furthermore this term 
is dominant in the ERBL region as we can infer from Fig. [6l 
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VI. CONCLUSIONS. 



In this paper we have defined the pion-photon vector and axial Transition Distribution Amphtudes using the Bethe- 
Salpeter ampHtude for the pion. In order to make numerical predictions we have used the Nambu- Jona Lasinio model. 
The Pauli-Villars regularization procedure is applied in order to preserve gauge invariance. 

We know from PCAC that the axial current couples to the pion. Therefore, in order to properly define the axial TDA 
(all the structure of the incoming hadron being included in A (x, S.,t)), we need to extract the pion pole contribution. 
In so doing, we found that the axial TDA had two different contributions, the first one related to a direct coupling 
of the axial current to a quark of the incoming pion and to a quark coupled to the outcoming photon and a second 
related to the non-resonant part of a quark-antiquark pair coupled with the quantum numbers of the pion. 

The use of a fully covariant and gauge invariant approach guaranties that we will recover all fundamental properties 
of the TDAs. In this way, we have the right support, a; G [—1, 1], and the sum rules and the polynomiality expansions 
are recovered. We want to stress that these three properties are not inputs, but results in our calculation. The value 
we found for the vector form factor in the NJL model is in agreement with the experimental result (2^ whereas the 
value found for the axial form factor is two times larger than in . This discrepancy is a common feature of quark 
models (T^]. Also the neutral pion vector form factor i^^^.^(i) is well described. These results allow us to assume 
that the NJL model gives a reasonable description of the physics of those processes at this energy regime. 

Turning our attention to the polynomial expansion of the TDAs, we have seen that, in the chiral limit, only 
the coefficients of even powers in ^ were non-null for the vector TDA. No constrain is obtained for the axial TDA. 
Nevertheless, the NJL model provides simple expressions for the coefficients of the polynomial expansions in the chiral 
limit, Eqs. ^ and 

We have obtained quite different shapes for the vector and axial TDAs. This is in part, at least for the DGLAP 
regions, imposed by the isopin relation ([55)1 . We have pointed out the importance of the non-resonant part of the qq 
interacting pair diagram for the axial TDA in the ERBL region. 

It is interesting to inquire about the domain of validity of the relations (j33p . These relations are obtained from the 
isospin trace calculation involved in the central diagram of Fig. 3. Due to the simplicity of the isospin wave function 
of the pion these results are more general than the NJL model and could be considered as a result of the diagrams 
under consideration. These diagrams are the simplest contribution of handbag type. 

The Transition Distribution Amplitudes proposed by the authors of [6] open the possibility of enlarging the present 
knowledge of hadron structure for they generalize the concept of GPDs for non-diagonal transitions. Calculated here, 
as a first step, for pion-to-photon transitions, these new observables should lead to interesting estimates of cross 
section for exclusive meson pair production in 7*7 scattering Q. 
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APPENDIX A: APPENDIX. 



In Section [IT] we have introduced the minimal kinematics definitions. The explicit expressions for the pion and 
photon momenta in terms of the light-front vectors are: 



2 2 



(1 - r + \^^^ + \ 



(Al) 
(A2) 



Here we have = [{p + p')/2f = ml/2~t/4, = {p' - pf , A-^'' = (0,Ai,A2,0) = (o,A-L,0^ and A^ = t. 

The polarization vector of the real photon, e, must satisfy the transverse condition, e.p' — 0, and an additional gauge 
fixing condition. When deriving Eq. ([5]), we need e.n = £+/P+ to kinematically become higher-twist, i.e. e.n 
when P"*" —^ 00. The standard gauge fixing conditions, = or e+ = 0, satisfy the previous requirement. In fact, 
all what we need is that the components of the polarization vector remain finite when P+ goes to infinity. 
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In Section IIIII we use the N JL model. We follow the notation of Ref . and we refer the reader to this paper for 
more details. The NJL model considers the lagrangian 



C = q{x){i^- Mo) q [x) + g [qqf + [qi^i^fqf 



(A3) 



where /xq is the current quark mass. As it is well known, the first consequence of the scalar interacting term is 
to provide a constituent quark mass, m, different from the current mass. Due to the point-like character of the 
interaction, this lagrangian is not renormalizable. We shall use the Pauli-Villars regularization in order to render the 
occurring integrals finite. This means that for integrals like the ones defined in Eq. (|14p . we make the replacement 



(A4) 



with m| = vn? + jA^, cq = C2 = 1, ci = — 2. Following Rcf. 8] the regularization parameters A and m are determined 
by fitting the pion decay constant and the quark condensate (in the chiral limit). With the conventional values 
[uu) = -(250 MeV)3 and = 93 MeV, we get m = 241 MeV and A = 859 MeV. 
The pion-quark coupling constant is given by 

(A5) 



with 



hip' 



(27r)^ (fe2 



2 



16^2 Z^ ^ 




le 



Am?- 

— 1 arctan 



Am'] 



(A6) 



1 



The numerical value of the pion quark coupling constant is g^^^ = 6.36 for the physical value of the mass of the tt"*" 
pion, and g^^^ = 6.71 in the chiral limit. 
In the NJL model, the vector form factor is 



it) 



SNr. 



9-nqq j / i\ 

T-^-^hip,? ) 



(A7) 



In order to calculate the form factors we need the expression for the three-propagator integral. In the particular case 
where p^ = rn^ and p'^ — 0, the expression for I^^pjp') is 



hip,?') 



d'^k 1 

(27r)4 (fc2 - m2 -I- ie) {{p + fc)2 - + ie) {{p' + fc)2 - m'^ + ie) 



2 



(4^)2 



dx 



:log 



x'^t + x{l - x)m1 



2m2 



VP 



(AS) 



y/p ^ x'^t + x{l — x)m1 — 2m'j + ^ 

with p = t^x"^ + x^{\ - x)'^m% + 2tx^{mlx{l - x) - 2m^), t = {p' -pf. 

The axial form factor involves the three-propagator integral as well, but also the two-propagator integral given by 

Eq.Ee] 



Ff (t) = 4A^, m g^gqV2 ( h ip, p' 



\h{^l)-h(t)\ 



(A9) 



The TT*^ — > 7*7 form factor can be obtained from the vector form factor through a isospin rotation: F^^*^ (t) 
\f2Fy^ (t) /m^. In the chiral limit and for ^ = all the considered form factors have the simple expression: 



F^ 



(0) 



V2 



F^ (0) 



Fl (0) 



V247r2/, 



o 2 



A2 



2A2 



0.192 X (1 - 0.108) GeV"^ = 0.171 GeV"^ 



(AlO) 
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The first coefficient of the right hand side is what is expected from the axial anomaly contribution to tt ^ 77 decay. 
The term between brackets has a small correction to the expected value of 1 due to the finiteness of the regularization 
masses. In the NJL model not only the quarks, but also the counter-terms run in the triangle diagram of the axial 
anomaly.In a proper renormalizable theory this correction disappears in the limit A — )■ 00. 
In the light-front calculations we need two kind of integrals. The first one is 



_(£k 

(27r)4 [{k -piY -m? + ie] [{k ~ ^2)^ - + ie 



1 1 ^ 
e(b - x)e{x - a) —- —- log 



m?{h ~ a)^ — {b — x){x — a){p2 — Pi)^ 
(47r)2 a — b ^ m'j{b — a)^ — (6 — x){x — a){p2 — Pi)^ 



(All) 



with b =max(l — p'^ /P^, 1 — P2 / P^) and a =min(l — p'l / , 1 — p^/P+). The second light-front integral needed 
in our calculations is the one defined by Eq. (jl4p . This integral can be performed by standard methods, obtaining 



{1-x)+C + VD 



\{l^x)+C- 



with 



D = 



- ^) - 2 (1 " ^) ) + "^,'(1 - C) (2m^e - t{l + i)) 



c= { 



\i\<x<\, 

-\^\<x<\i\ ^>o 



-[i^x) 



.22C(1"C) 



-Iel<a:<|e| 



As an aplication of the expression (jAlip we can evaluate the PDA. From its definition in Eq. 



4 TVe m g^q 



d'^k 

(27r)'* [fc^ — m? + ie] [{k — p)^ — w? -\- it 



This integral is of the type of (jAlip with a = and 6=1. Therefore 



(A12) 



we have 



(A13) 



{x)e{i~x)Y<^,\og 



TO J (1 — x) X 
TO^ [1 — x)x 



(A14) 



As it is expected for the PDA, x runs from to 1. We can test our result verifying that in the chiral limit (j){x) = 
6 (x) (1 — x) , as it is well known for the NJL model. This peculiarity of the NJL model in this limit is present also 
in the quark valence distribution, which becomes as simple as q{x) = 6 (x) 6{1 — x). In Ref. [l3] is discussed how 
this is, in fact, a quite reasonable approximation of realistic models. 
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